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ABSTRACT 


The Volume dependence of the Griineisen parameter has been used to 
calculate the pressure dependence of melting temperature for MgO and hep 
iron using Lindemann law. The Volume dependence of the Griineisen 
Parameter has been determined using gamma volume y(V] relationship 
due to Al'tshuler and reciprocal gamma-volume y(V] relationship due to 
Srivastava et al. The calculations are performed between pressure range 
55GPa-330GPa. The values of Melting temperature are calculated at 
different pressures and compared with the available data reported in the 
literature. 
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2. VOLUME 
PARAMETER 

Anderson et al.[l] have investigated an empirical 
relationship between y and V for hep iron using the 
laboratory data. This relationship is writer as follows 

y = 1.05051nV — 0.2799 (2] 

Anderson et al. [1] reported Eq. (2] as a new result based 
on the experimental data for hep iron, and used it in the 
Gilvarry law to determine values of T m at different 
pressures. It should be emphasized here that Eq. (2) is not 
physically acceptable. 

Dorogokupets and Oganov [15] have used the Al'tshuler 
relationship [6] for y(V] in order to determine the 
thermoelastic properties of solids at high pressures and 
high temperatures. According to this relationship, we can 
write 

Y = Yoo + Oo ~ Xoo) Qr) (3) 

Where n is a constant depending on the material. y 0 and 
y^are the values of y at zero pressure and in the limit of 
infinite pressure respectively. V 0 is the volume V at zero 
pressure. On differentiation Eq.(2] gives the second order 
Gruneisen parameter 

• -iSM 1 -?)* t4) 


1. INTRODUCTION 

Anderson et al. [1] have presented an analysis of the high- 
pressure melting temperature of hep iron using three 
different approaches viz . the Lindemann -Gilvarry law 
[2,3] taking in to account the volume dependence of 
Gruneisen parameter y the Clausius-Clapeyron equation 
accounting for the liquid state of melting , and the Poirier 
dislocation-mediated melting [4,5]. The Lindemann- 
Gilvarry law, most frequently used in the recent literature, 
is writen as follows 


d\nT m 
d In V 




(i) 


In order to determine values of melting temperature T m at 
different compressions or pressures by integrating 
equation (1], we need to know y as a function of volume V. 
there have been various attempts [1,6-7] to develop 
formulations for y(V). There have been various attempts 
[6-11] to formulate expressions for y(V] as well as y(P]. 
The most critical examination of these formulations can be 
made by applying the thermodynamic constraints for the 
higher-order Gruneisen parameters [12, 13-16] at 
extreme compression in the limit of infinite pressure. In 
the presented study we have selected MgO which is an 
important ceramic material and geophysical mineral [1], 
and hep iron which is the main constituent of the earth 
core [8,12,16]. 
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And further differentiation yields the following expression 
for the third order Griineisen parameter 

! _ d\nq _ Yoon ^ 


d In 


Eq. (4) reveals that q tends to zero when y tends to y^. 
This is consistent with the thermodynamic constraint 
[2,3,13]. In the limit of infinite pressure, Eq. (5] gives 

( 6 ) 


2-00= n 


And therefore 
Ay =A 00 7 oo= constant 


(7) 


It has been found [2,3,13-16] that Aandy both decrease 
with the increase in pressure. However, Eq.(7] which is 
based on Eq. (2) is in contradiction to this finding 
according to Eq. (7], A must increase with P when y 
decreases. Thus it comes evident that Eq. (2) is not valid. 

In order to rectify this shortcoming, attempts have been 
made to use the following expression for reciprocal 
gamma [17-11] 


i _ j_ /jl _i_\ f^_\ n 

Y Yoo Vy 0 YooJ \V 0 J 


Yoo \Yo YooJ \V 0 J 

Eq. (8) on successive differentiation yields 


dlny (y \ 

q = -= n-1 

^ d\nV Vyoo / 


and 
A = 


d\nq_ny 
d\nV yoo 


( 8 ) 


(9) 


( 10 ) 


Eq. (9] is consistent with the constraints that q tends to 
zero when y tends to y& Eq. (10] gives 2^ = n and 


- = —= constant 

Y Yoo 


( 11 ] 


Eq.(ll] is consistent with the finding that A and y both 
decrease with the increase in pressure. For determining 
values of the third order Grniieisen parameter A at limit of 
infinite pressures in case of the earth lower mantle and 


core, Stacey and Davis [2] scaled — to — .This is an 


A 

Aoo yoo 

agreement with Eq. (11]. A more critical test of an 
expression for y (E] can be provide using the definition of 
the fourth order Griineisen parameter ^ given below [16] 




din A 
din 


( 12 ) 


On integrating Eq. (13],we get 


And further integration gives 


Y o 


A 

Aoo 


e-ooT 


gO 


(15) 


(16) 


It has been found [18,11] that if we assume 

2o — 2 oo = <7 0 (17] 

Then Eq. (16] reduces to Eq. (8], and Eq. (14] becomes 
compatible with Eq. (9] to (11]. Eq. (17] is not satisfied by 
the available data [2,16]. Eq. (14] on differentiation yields 


(18] 


In the present study we use Eq. (14], (15], (16] and (18] to 
obtain y, q, A and ^ for MgO and hep iron. 

3. LINDEMANN LAW OF MELTING 

Lindemann law is represented by Eq. (1]. We prefer Eq. 
(8] over Eq. (14] to be used in Eq. (1] because Eq. (8] can 
be integrated analytically to obtain values of T m at 
different compressions. It is difficult to integrate Eq. (14] 
using analytical methods. Eq. (1] and (8] taken together 
yield on integration 


2 yoo 


Tm _ (Y_\ Aoo 7M 

T m o W vv 


-2yoo+- 


(19) 


We have used eq. (19] to determine values of T m (V). The 
results are then obtained for 7 m (P) using the pressure- 
volume relationship for hep iron[9] based on the seismic 
data [8]. 

4. RESULTS AND DISCUSSIONS 

For hep iron, we take y 0 = 1.83, y^ = 1.33, K 0 = 5.0 ,q 0 = 
1.18 from Stacey and Davis [8]. In Eq. (13], T mo , V 0 and y 0 
correspond to the reference or the starting point. For hep 
iron the starting point is at about 55GPa and T mo = 
T m (55GPa) = 2790 K [1]. Data have been used to obtain 
V(P]/V(55GPa] with the help of the Stacey reciprocal K- 
primed E0S[8,9]. 

We use the Stacey K-primed equation of state, which is 
based on Eq. 


(-) 


i 

Koo 


( 20 ] 


Eq. (11] and (12] yield 

din A _ dlny 
d\nV d\nV 


(13) 


According to Eq. (13] the fourth order Griineisen 
parameter ^ is equal to the second order Griineisen 
parameter q. This result based on Eq. (8] is not consistent 
with the results for ^ and q recently reported by Shanker 
et al. [16]. It should be emphasized that Eq. (8] is a 
simplified version [18, 11] of the Stacey- Davis 

formulation [2] for y(E]- According to the Stacey- Davis 
model we can write 


A — Aw + (2 0 


go 


(14) 


And can be written as follows; 

, p 


± = ±+( 1 -?k) 

K' K n V K n J 


K 


On integration, Eq. (21] gives 

P \-K 0 /Koo 


K 

«0 


= 


Further integration of Eq. (23] gives 


( 21 ) 


( 22 ) 


(23) 
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Table (1] values of volume compression V(P)/V(55GPa) at different pressures for hep iron based on the seismic data and 
the Stacey equation of state [8,9], Griineisen parameter y calculated from Eq. (3) at n=3.66 and from Eq. (8) at n=3.14 and 
2 . 66 . 


T„, 

P(GPa) 

V(P)/V(55GPa) 

Y Eq.(3) 

YEq.(7) 

2790 

55.0 

1 

1.571 

1.557 

3014 

68.8 

0.9686 

1.544 

1.532 

3195 

80.5 

0.9452 

1.526 

1.515 

3476 

100.0 

0.9117 

1.502 

1.493 

3892 

131.1 

0.8679 

1.473 

1.467 

4196 

155.6 

0.8394 

1.457 

1.452 

4551 

186.0 

0.8092 

1.441 

1.438 

4971 

224.2 

0.7774 

1.426 

1.424 

5240 

250.0 

0.7587 

1.418 

1.417 

5736 

300.0 

0.7275 

1.405 

1.406 



5. CONCLUSION 

The results for T m at different pressures and compressions 
are given Table. Values of T m for Mgo and hep iron 
determined from Eqs. (14] and [15] are in good agreement 
with those obtained from the Stacey -Irvin model [8] 
based on the fundamental thermodynamics [1]. The most 
important conclusion is that the Lindemann-Gilvarry law 
of melting is valid for hep iron and consistent with the 
thermodynamic constraints for the volume dependence of 
the Griineisen parameter at extreme compression. At 330 
GPa, the pressure for the Earth inner core-outer core 
boundary, the melting temperature for hep iron is very 
close to 6000K. The importance of infinite pressure or 
extreme compression behavior of materials has been 
discussed in details by Stacey and Davis [8] pointing out 
the usefulness of infinite pressure constraints in 
determining the properties at finite pressures. Values of y, 
q, X and ^ at and hep iron calculated from the two 
formulations compare with each other. 
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